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ABSTRACT 


The analysis and control of errors in precision 
mechanisms is extremely important in order to obtain the 
desired level of accuracy. In this work three types of 
mechanisms, namely, Geneva mechanisms, cam follower system 
and gear train are considered for mechanical error analysis 
and synthesis. Both deterministic and probabilistic pro- 
cedures are presented for predicting the effect of various 
tolerances and clearances on the acceleration and jerk of 
external and internal Geneva mechanisms. The effect of 
changes in geometrical and other system parameters on the 
kinematic and dynamic response of a cam follower system is 
found stochastically. The probability characteristics of 
the integrated gear train position error of gear train are 
determined using the method presented by Michalec. The 
problem of allocation of tolerances and clearances (3ar 
errors) in mechanisms is formulated as a constrained optimi- 
sation problem and nonlinear programming techniques are used 
for the solution. This procedure distributes the tolerances 
so as to minimize a measure of the manufacturing cost while 
satisfying the specified requirements on the deviation of 
the output parameter of the mechanism. The present study 
is expected to offer an unified method of error analysis 
and synthesis for all mechanical systems. 



CHAPTER 1 


INTRODUCTION 

1 . 1 IMPORTANCE OF EEROR ANALYSIS 

The development of science and technology is 
closely connected with the creation of modem high preci- 
sion mechanisms, machines and other sensitive devices 
which produce the parameters of motion with high accuracy, 
Mechanisms are applied widely in measuring devices, compu- 
ting devices, recording and registration equipment, optical 
instruments, automatic-controlling devices and similar 
others. 

It is well known that in a real mechanism, the 
actual dimensions will he different from the theoretical 
ones. The difference will he caused hy the limitations of 
the manufacturing processes, poor workmanship, inaccuracies 
of mechanism assembly and the clearances in the hearings 
of kinematic pairs which provide relative displacements of 
links. If the tolerances and clearances exceed a particular 
limit, then the very purpose of the mechanism will he lost. 
For example, if the tolerances and clearances exceed a 
particular limit in a Geneva mechanism, then the roller 
may enter the slot with a large value of jerk, which will 
ultimately result in poor performance of the mechanism. 

In the case of internal combustion engines with valves, 
the lift of the valves, which is governed hy a cam-follower 
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system, should "be known very accurately so as to clearly 
define the closing and opening of the valves. Otherwise 

it may result in the inefficient use of the fuel. Similarly, 

/ 

in the case of gear-train used in precision mechanisms, the 
various mechanical errors have to he controlled in a proper 
manner in order to obtain the desired level of precision in 
terms of motion, power transmission and dynamic response. 

1.2 LITERATURE SURVEY 

Several attempts have been made to predict and 
control mechanical errors in mechanisms in the past. 
Hartenberg and Denavit \ 1 j presented a deterministic method 
for analysing the error in four-bar planar mechanisms. 

This method does not take into account the effect of clear- 
ances. A deterministic method has been presented by 
Kolhatkar and Yajnik \ 2 \ for estimating the effect of 
clearances on the output of a mechanism. Tuttle |3l used 
vectorial addition to analyze errors in cam operated devices 
and four-bar linkages. This method gave the approximation 
only of the error in a non-continuous system, and is not 
applicable to error analysis of cam-follower mechanisms. 
Rothbart jM used a finite difference method to obtain 
the effect of surface roughness on follower acceleration 
(deterministically). The effect of manufacturing errors 
on the integrated gear train position error of a gear 
train was considered deterministically by MLchalec j 5 \ • 
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In this work the total error was computed by considering 
the contributions of individual errors as additive. Hence 
this method gives a conservative estimate of the total 
error. 

A statistical method for the mechanical error 
analysis of planar four-bar mechanisms was offered by 
Garett and Hall 1 6 j . This method yielded some useful 
results and paved the way for further research in that 
direction. This method was computer oriented and a lack 
of rigorous mathematical treatment makes it unsuitable 
for use in synthesis. Dhande and Chakraborty \ 7 \ have 
provided a stochastic model of the four bar linkage con- 
sidering tolerances and clearances. The mechanical error 
was analyzed for the three-sigma band of confidence level. 
The problem of allocation of tolerances and clearances for 
a four-bar mechanism in an optimal manner was also consi- 
dered by them. Rao and Ambekar j8j presented the error 
analysis of spherical function generating mechanisms by 
using a probabilistic approach. Hiande and Chakraborty 
employed a stochastic model by considering the kinematics 
of spatial cam mechanism, to analyse the three-sigma band 
of error in the follower displacement showing maximum and 
minimum errors |9j . Kim and Newcombe jlO| extended the 
method of reference f9j to enable an analysis of the error 
in velocity and acceleration due to any given manufacturing 
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tolerances on the cam-profile. A probabilistic error 
analysis of gear trains was presented by Michalec j $\ . 

By idealizing the gear train as a weakest link kinematic 
chain, Bao presented a method of designing gear trains 
with respect to bending strength and surface wear resis- 
tance jllf . 

A model, which has been named as an impact pair, 
has been made by Dubowsky and Freudenstein |12j to analyze 
deterministically, the dynamic response of mechanical 
system with clearances, by considering the displacement- 
forced motion, which is quasi-harmonic. While discussing 
this paper (at the end of Ref. |l2( ) Erdman and San dor 
commented that at high speeds, the varying inertia forces 
in the mechanism as well as the strains due to elasticity 
of links, will certainly affect the forces on the impact 
pair. In other words, there will not only be displacement 
forced motion, but other forces, which will be functions 
of parameters such as length and stiffness of links, confi 
guration, distribution of mass and input speed of the 
mechanism, will also be there. The idea of considering 
the effect of errors on the dynamic response of cam- 
follower systems, presented in Chapter 3 of this thesis, 
germinated from these comments. 
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1.3 SCOPE OF THE PRESENT WORK 

In this work, the analysis and synthesis of 
mechanical error is considered for the following mechanisms : 

(1) External and internal Geneva mechanisms 

(2) Cam-follower systems 

( 3 ) Gear-trains 

In the case of Geneva mechanisms, firstly, the 
effects of tolerances on acceleration and jerk have been 
studied deterministically. After that, taking into account 
the clearances as well, a stochastic approach has been used 
and the three-sigma band widths of acceleration and jerk 
are plotted. Finally, the tolerances and clearances are 
allocated optimally with upper bounds on the standard 
deviations of acceleration and jerk. 

In the investigation of cam-follower system, a 
probabilistic approach has been employed to analyze the 
effect of manufacturing errors in cam profiles and other 
geometrical parameters on the kinematic response. The 
variations in other system parameters such as equivalent 
mass, stiffness and damping are also considered to find 
the probabilistic characteristics of the dynamic response 
of the cam-follower system. The synthesis problems, invol- 
ving the optimal allocation of errors, have also been 
solved by placing upper bounds on the standard deviations 
of the responses. 
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In the case of gear-train problem, the probabilis- 
tic error analysis presented in Ref. j5j has been used 
while solving the synthesis problem. 


1 .4 METHODOLOGY 

1,4.1 For Deterministic Error Analysis: 

For any mechanism, the equation relating the 
system parameters (a-| , a 2 , ..., a n ), input variable (9) 
and the output variable (0) can be expressed in a general 
form as 

f (a.,, a 2 , a^j 9 ; 0) = 0 (1.1) 

If Aa i (i = 1 , 2, . . . , a n ) represents the 
error in the system parameter a^, the error in the output 
variable ( A 0) for any fixed value of 9 can be expressed 
as 

f = f ( a^ + a^ , a 2 + 0A a 2 j • • • + -A 5 9 , 0 ■+■ /A 0 ) = 0 

(1 . 2 ) 

If a ^ represent the tolerances or clearances, 
then A 0 is called the mechanical error. Equation (1.2) 
can be expressed using Taylor' s series expansion as 


f = f + df + ♦ . . 


f ( a-j , a 2 , . # • , a -Q 5 9 5 0 ) 


•_^ n 'bf 

•^=1 SR 


na ± + y|l 


C 3^ ^ 3^ j * p • 9 9 ) 


f f ? 0) 


0 4- *• • 


0 


(1.3) 
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For small values of A a i> df = 0 and hence 
the mechanical error ( -A.0) in the mechanism can be found 
as 



where the notation j has been used to denote the evalua- 
tion of the quantity at (a^, a 2 , . .., a R ; ©). 


1.4.2 For Probabilistic Error Analysis: 

The various tolerances and clearances in the 
mechanism have been assumed to be random variables. In the 
absence of any information, these random variables can be 
taken to be normally distributed. As the output variable 
(0) of the mechanism (like displacement, velocity and 
acceleration) is a non-linear function of the system para- 
meters (including those for which tolerances and clearances 
are specified), it can be expressed as 


0 — 0 ( a ^ , a 2 j » • « , a^ | 0) ( 1 • J? ) 

Thus 0 will also be a random variable. As the 
thecory of random variables is well established for linear 
functions, the Taylor’s series expansion about the mean 
values of the random parameters is used to linearize the 
function 0. Since the sum of several random variables 
following normal distribution is also a normal variate, ; 

0 can also be treated as a normally distributed variable, 



8 
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whose moan ( /-^j) and standard deviation ( CT^) are 
given by 


J\ & ( > Ma 1 » > P a 2 » *** 5 / A a n 


e) 


( 1 . 6 ) 


and 


4-\i < a *> 


i=1 


( 


■ > / A a 1 » ^a 2 » 


a. 


e ) 1 


2 

*7~a, 


where 


^ ^ aj 


1 

n J (1.7) 
denote the mean and the standard devia- 


l “i 

tion of a^, respectively. Although the output variable 0 


has been assumed to be normally distributed on the basis of 
linearization, it can also be justified from central limit 
theorem M . 




1.4- .3 For Synthesis Problem: 


The synthesis problem, involving the allocation 
of tolerances and clearances in the mechanisms, has been 
solved as a nonlinear optimization problem. As the manu- 
facturing costs increase very rapidly as the tolerances and 
clearances to be maintained, approach zero, the objective 
function for minimization can be taken as the sum of the 
reciprocals of the tolerances and clearances. This func- 
tion will represent, though not the exact manufacturing 
costs, at least a measure of the manufacturing costs. 

By taking the tolerances and clearances as the 
design variables x^, upper and lower bounds are placed 


// 




on each x^ depending on the manufacturing limitations. 



Further, bounds are also placed on the deviation of the 
output variable. Thus the aim of the optimization problem 
becomes one of minimization of the cost, that is, maximiza- 
tion of the tolerances and clearances, keeping the devia- 
tions of the output variables well within the specified 
limits. It is important to note that the results of optimi- 
zation can not be applied directly for practical designs^ v ’ 
The tolerances and clearances given by the optimization 
method have to be rounded off to nearest feasible discrete 
values . 

For the solution of the constrained optimization 
problem, the interior penalty function has been chosen, 
coupled with variable metric method^of unconstrained mini- 
mization and the cubic interpolation method of one dimen- 
sional minimization |l4-f . If the optimization problem is 

stated as, 

Find X, which minimizes f ( X ) 

, ( 1 . 8 ) 
and satisfies g. ( X ) 4o, j = 1, 2,' 3,,.., m 

J 

the following iterative process can be used in 
the interior penalty function method. 

(i) Start with an initial feasible point X^ satisfying 
all the constraints with strict inequality sign, 
that is, gj (X 1 ) for j = 1, 2, ..., m, 

and an initial value of r^ 0. Set k = 1. 
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(ii) 


(iii) 


(iv) 


(v) 


m 

Minimize 0 ( X, r v ) = f ( IT ) - r, ^ — , 

k k ^i gj(i) 

by using the variable metric method and obtain the 


solu 'cion 


4 


Test whether X^ is the optimum solution of the 
original problem. The convergence criterion is as 
follows : 

(a) The process is assumed to have converged whenever 
the relative difference between the values of the 
objective function obtained at the end of any two 
consecutive -unconstrained minimizations falls 
below 0.00001, or 

(b) The process is to be terminated when the diffe- 

— . 

rence between the optimum points X^. and 
becomes very small as given below 

IcAltXjJ o.ooooi ( 1 . 9 ) 

where AX = X^ - X^_^ and ( -A X is the 
ith component of the vector AX. 


Find the value of the next penalty parameter, r k+ ^ r 


as 


^+1 = 0>1 ( V 

Set the new value of k = k + 1 , take the new 
starting point as X^ = X.Z and go to step (ii) 


(1 . 10 ) 
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1 .5 ORGANISATION OF THE THESIS 

The analysis and synthesis of external and inter- 
nal Geneva mechanisms is given in Chapter 2. Chapter 3 
deals with the analysis and synthesis of the kinematic as 
well as the dynamic response of the cam-follower system. 

The synthesis of the gear-train is investigated in Chapter 
Conclusions and recommendations have been included as 
Chapter 5* Appendix A contains the partial derivatives 
of acceleration and jerk of a Geneva mechanism, needed 
for the synthesis problem of Chapter 2. The probabilistic 
analysis of gear train errors constitutes appendix B. 
References are given at the end of the thesis. 



CHAPTER 2 


GENEVA MECHANISMS 

2.1 GENERAL DESCRIPTION: 

Geneva mechanisms have long "been popular as a 
means for producing positive incremental motion because of 
their simplicity, both in design and construction, which 
makes them relatively low cost indexing devices. In a 
Geneva mechanism, the input shaft rotates at a constant 
speed while the output, shaft turns through a suitable angle 
and then dwells for a given interval before repeating an 
identical motion. The output shaft always rotates in the 
same direction. 

(i) External Geneva Mechanism: 

Figure 2.1. a shows an external Geneva mechanism 

where 

D = driving or input wheel, 

F = follower or output wheel, 

L = locking arm, and 
2 = indexing period of the driver. 

The centre line of the slot should be tangential 
to the pitch circle of the pin at the points of engagement 
and disengagement. The minimum number of slots is three 
and in a practical mechanism it is advisable to make the 
number of slots higher than the absolutely minimum number. 



13 


More number of slots give less noise il 7 \„ Unfortunately, 
the locking of the wheel during the dwell period becomes 
less secure with larger number of slots. The mean accelera- 
tion of the follower (that is when no tolerances and clear- 
ances are considered) is governed by the number of slots. 

The indexing period of the Geneva mechanism shown in 
Fig. 2.1. a with four slots is equal to 90°. 

(ii) Internal Geneva Mechanism: 

In the internal Geneva mechanism shown in Fig. 
2..1.b, D = driving wheel, P = pin or roller on the driving 
member, F = follower or output wheel, L = locking arm, and 
2 u/ = indexing period of the driver. The indexing period 
of this internal Geneva mechanism is 270°. 

External Geneva drives have short indexing periods 
and hence the accelerations will be large because the 
acceleration is inversely proportional to the square of the 
indexing period. An internal Geneva drive has a long 
indexing period and hence the magnitude of the dynamic force 
(which is proportional to acceleration) is considerably 
lower than in the case of an external Geneva drive having 
the same number of slots. The character of the accelera- 
tion curve for an internal Geneva mechanism limits . its 
application to low speeds only. The values of jerk are 
much smaller in the case of internal Geneval mechanism 
than in the case of external Geneva mechanism. The motion 
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of an internal Geneva mechanism near the central position 
can he considered to he uniform for all practical purposes 
as shown in Fig. 2.1.c il8 \ . 

2.2 DETERMINISTIC ERROR ANALYSIS OF EXTERNAL 
GENEVA MECHANISM: 

The following nomenclature is common to hoth 
external and Internal Geneva mechanisms : 

a Q = centre distance between the driving wheel and 

the follower with no tolerances and clearances, 
b Q = distance between the centre of the pin and the 

driving wheel with no tolerances and clearances, 
a = nominal centre distance = a Q , 

b = nominal distance between the centre of the pin 

and the driving wheel = b c , 
c = width of the slot 

/\a , &b = tolerances in the link lengths a and b 
respectively, 

A c = clearance in the slot width 
b + fib = b^ = actual length between the pin and the 

driving wheel, 

a + A a = a_ 1 = actual length between the pin and the 

follower 

(X = input angle 
P = output angle 
Ac = c ^ 
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% - 


A = 


ACC = 


JEM 


N 


p+ 

1 


d 2 p 


d<X 

,2 


2 ~ - angular acceleration 


dt 


( 


2^ N \ 2 


60 


) . A 


J = 


d 3 g > 1 


c?^ 

l3j 

i 

dt- 


d3g 1 

.3 


2 A N 


) . J 


^ 60 

speed of the driver (r.p.m.). 


Figure 2. 2. a shows, schematically, the roller 
in an intermediate position in the slot. Here it is assumed 
that the roller always remains in contact with the slot and 
that it never floats in it. From the geometry of the mecha- 
nism, the following relation can he derived 2 


Sin P 


b sin-x 

o • 

° (a o + b o - 2 a o b o Cos ^ o )1/2 


( 2 . 1 ) 


Figure 2.2.b shows the mechanism with tolerances 
and clearance. The following geometrical relations can be 
derived 


Sin (<X + ft + A f’b ) = A C/x 


( 2 . 2 ) 
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x ± b 4- A b 
Sin ( + A fl ) 


a + A a 

Sin 7T - (On + f3+ A G> % 


(2.3) 


Ea. 


Substituting the value of x from Eq. (2.2) into 
(2.3) and solving, it can be shown that 


. A 

Sin •" 


b^ Sin^ 

(a 2 + b 2 - 2 a 1 b 1 Coso< ) 1/2 


where 


'1 


(a 1 - b^ Cos«X ) 

(a 2 + b^ - 2 a 1 ^ Cos 'X ) 


= a + A a 
= b +A b 
= A c 


(2.4) 


The second and the third derivatives of the output angle, 
with respect to the input angle can be derived as 


d 2 |3 1 


B. 


! i.2 

(A 1 + ) (A 1 - B 1 ) 


d* 


1 1 ~ (A^ + B 


PP 75 + ft - U, ~ ip 2 ! 372 

(2.5) 


where 


A 1 

a! 


B 


B, 


b 1 Sinc^/T 


b^ Cos c</T - a n b 2 Sin 2 o< /T 3 
Cl (a 1 - b 1 CosO< ) / T 2 


b 1 c 1 SinCX/T £ 


2 a^ b^ c^ Sin ok 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


(a^ - b^ Cos c* ) / T 


(2.9) 
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B . 


T 


a 3 A | 

do(3 


Here 


A . 


IM 


B„ 


t » 


~ SinCK. /T - 3 B 2 Sine* Coso< /T 3 
+ 3 B 3 Sin 3 <x /T^ (2.10) 

B 1 c 1 Cos <x /T 2 - [ 4- a 1 B 2 c 1 Sin 2 ^ •; 


1 


+ 2 a.j B^ Coso< (a^ - B^ Cos t* ) j| / 

T* + 8 a 2 B 2 Sin 2 o< (a 1 ~ B 1 Coso^ ) / T 6 


2 ? 1/2 
(a^ + Td^ - 2 a 1 B 1 Cos ) 


( 2 . 11 ) 

( 2 . 12 ) 


tit It! 


It It. 


(A 1 + B 1 ) / Y + 3 (A 1 + B 1 ) (A 1 + B 1 ) (A 1 + B 1 )/Y 3 

+ (A* + B^j ) 3 / Y 3 + 3 (A 1 + B 1 ) 2 (A^ + b' ) 3 / Y^ 

(2.13) 

r 2-, 1/2 

[ 1 - (A 1 + B 1 ) J (2.14) 

Cosc^/T + (4 a^ B 2 Sin%* - 3 a 1 B 2 Cos 2 c<)/T 3 

+ 18 a 2 B 3 Sin 2 «X Cosot./ - 15 a 3 B^SlnVl? 

(2.15) 

- B 1 Sine</ T 2 + 2 a 2 b^ 3inc</ 

, p Cos<* 4 p p 

- 14 a^ B^ c r sin <=</ T + 24 a^ B^ sin 0C • 

CoscX (a^ - B^ Cos <X ) / + 24 a 2 B 3 Sin^*/ 


- 48 a 3 B^ Sin 3 ex . (a^ - b^ Cos 04, ) / T 1 


,8 


( 2 . 16 ) 
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The individual effect of tolerances and clearance on the 
acceleration and the jerk is shown in Figs. 2.3 to 2.8 for 
an external Geneva mechanism having 4- slots with a crank 
length of 10 cm. From Figs. 2.3 to 2.5, it can be observed 
that the effect of A c on the acceleration is much less than 
the effects of A a and A b . Hence a greater variation can 
be tolerated in the slot width compared to the variations 
in the crank length and the centre distance. It is to be 
noted that in these figures, the points lying on the Y - axis 
are singular because the Y - coordinate (which represents 
the percentage change in acceleration) blows up for zero 
mean acceleration (denominator) ato<= 0. Figures 2.6 to 
2.8 show that the variation of jerk does not follow a 
definite pattern. Here also the tolerances in the crank 
length and the centre distance can be seen to be more impor- 
tant. The maximum percentage changes in acceleration and 
jerk for various values of A a, A b and A c are given in 
Table 2.1. 


2.3 DETERMINISTIC ERROR ANALYSIS OF INTERNAL 
GENEVA MECHANISM: 

Figures 2.2. c and 2.2.d show an internal Geneva 
mechanism without and with tolerances. By proceeding as 
in the case of external Geneva mechanism, it can be 
derived that 
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(2.17) 



Sin ^ 


ih 

4 *? 


where 


b Q Sin os 

(a 2 + b 2 + 2 a Q b Q Coso0 1//2 
b 1 Sin os 

(a 2 + b 2 + 2 a^ b,j Cos c*. 
c H (a., + b. Cos cX ) 

+ - V— U r — 1 (2.18) 

(a^ + b^ + 2 a 1 b^ Cos cx ) 

2 

A^ + B^ (A^ + B^)(A^ - B^ ) 

F^TT 75 + n^^vT 75 

(2.19) 


= b 1 Sin cx/ T (2.20) 

B 1 * c 1 (a 1 + b 1 Cos c*i ) / t 2 (2,21) 

A* = b 1 Cos CX/ T + &1 b 2 Sin 2 cx. / T 3 (2.22) 

A^ =~b 1 Sin cx / T + 3 b 2 Sin ox Cos^K / 

+ 3 a 2 b 3 Sin 3 cX / (2.23) 

B 1 = c 1 (a 1 + b 1 cos cX ) / T 2 (2.24) 

B^j = — c 1 b 1 Sin CX/ T 2 + 2 a 1 b^ c 1 Sine*.. 

(a., + b 1 Cos<x ) / T h (2.25) 


b” = - b 1 c 1 Cos OS/ T 2 - 4 a 1 b 2 c 1 Sin^-X/ T 1 * 

+ 8 a 2 b 2 c 1 Sin 2 cx(a 1 + b 1 Cos <X ) / T 6 

+ 2 a 1 b 1 c 1 CosO< (a^ + b^ CosoO/T* 4 " ( 2 . 26 ) 

1/2 

T = (a 2 + b 2 + 2 a 1 b., Cosoc) (2.27) 
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do , 3 


where 


Y = 


A. 


tu 


E, 


IT! 


= (A*" + b”' ) / X + 3 (A 1 + B 1 ) (a' + b' )(a’' + b")/X 3 

3 2 2 

+ (a* + b* ) ' / y 3 + 3 (A 1 + b 1 ) (a’ + B* ) / y^ 

( 2 . 28 ) 

2 - 11/2 

j 1 - (A 1 + B 1 ) / (2.29) 

- b 1 Cos ex. / T + (3 a 1 b^ Cos^i - 4 a 1 bi~ Sin§x. ) 
/ T 3 + 18 b 3 SinlcCosCX/ T^ 

+ 15 a 3 hj Sintb< / T 7 (2.30) 

b 1 c 1 Sincx,/ T - 2 c 1 b 1 Sin o(/ T* 

- 14 a 1 c 1 b^ Sin <x . Cos <X / T* 

+ 24 a^ b^ Sine**. Cos d^a^ + b^ Cosoc )/ 

- 2b b 3 Sin 3 

+ 4-8 a 3 b 3 (a 1 + b-j Cos c< ) / T^ (2.31) 


Numerical results have been obtained for the changes in 
acceleration and jerk due to variations in a, b and c for 
an internal Geneva mechanism having 4- slots with a crank 
length of 10 cm. The results are given in Tables 2.2 to 
2.f.. A summary of the maximum percent changes in accelera- 
tion and jerk obtainable for various values of ^ a, A b 
and A c is given in Table 2.6. It can be observed that 
contrary to the results of the external Geneva mechanism, 
here the effect of Ac (that Is the clearance in the slot 
width) on the acceleration is greater than the effects of 
A a and A b. 



21 


2.k PROBABILISTIC ERROR ANALYSIS CONSIDERING 
TOLERANCES AND CLEARANCES FOR THE 
EXTERNAL GENEVA MECHANISM 


Figure 2.9 shows a probabilistic model of the { 
external Geneva mechanism by considering both tolerances 
and clearances as random variables. Some clearance between 
the pin of one link and the race of the adjacent link is 
always present in practice. Here r^ 2 is the radial clear- 
ance between the output wheel (follower) and its shaft and 
r^-j is the radial clearance between the input wheel (crank) 
and its shaft. Obviously, the clearance zone between the 
roller and the slot is rectangular. r^ shows the width 
of this rectangular clearance zone. A rectangular coordi- 
nate system (x. ., y. .) is so chosen that the x . . axis is 
coincident with the centerline of the link ending in the 
race. For most of the practical problems it can be safely 
assumed that the probability of the pin axis lying at any 
point is same for all the points inside the clearance zones. 
The equivalent linkage is shown by dotted lines in Fig. 2.9. 
From the geometry of this figure, it can be seen that 


a e = (a i + x i2 )2 + y ?2 

2 

bf = (b 1 + * 31 ) + yf, 

°e = ^23 


( 2 . 32 ) 


Since (a 1 +x 12 )>»y 12 and (b 1 + x^ ) , 

hence a Q o*, + x 12 and b g b 1 + x^-j . In the pre- 

vious section, the expressions for acceleration and 
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jerk have been derived as 


Acceleration 

= A 

(&-] > 

V 

c 1 ,0< ) 

Jerk 

= J 

( a -| ) 

V 

c -j ) 


(2.33) 


By considering the equivalent linkage with tolerances and 
clearances, the new expressions for acceleration and jerk 
can be written as 


Acceleration = A(a,b, c,(X) 

© G @ 


Jerk 


= J ( a e > c e’ ) 


(2.3^) 


X * © # j 

Acceleration = A (a +Aa + x^ 2 , "b +Ab + x^ , J^yOc) 
Jerk = J (a +Aa + x 12 , ^ +J ^ 1 ° + » y 23 ,<:> ^ 

(2.35) 

In the above expressions, A a, Ab, x.j 2 , y 23 and are 
random variables. 

If y"j_j are the coordinates of the pin axis, 

the probability density function is given by 


f (x 12 , y 12 ) = 


r JL 

r 


V 2 
r 12 


0 

1 


f U 31 , y 31 ) = 


_ 2 
7\r 3 i 


0 


,2 

"12 


x. 


+ 4 > 


2 

•v» 

r 12 

+ jf 2 

> 

r 2 

r 12 

2 


2 

+ y 31 

r 31 

+ y 31 

> 

r 2 

r 31 


( 2 . 36 ) 


(2.37) 
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The upper bound on the length of the rectangular 
clearance zone is (2 r^ + 2 r^). Hence 

i r 


t (■ X 23> ^23^ 


2(r 12 + r “ J 2. » 


0 


\ (~ 2 r 12^ x 23^- 2 r 31^ ^ 
( - r 23^ y 23^ r 23 } 


(2.36 


otherwise 


K 


If X 1 and are two random variables, then 

mom 

rr 


the correlation moment of X^ and X 2 is defined as 


Z 1 X 2 


Hence 


J 


x 


1 x 2 f (x , x ) dx dx 2 


(2.39) 


K. 


X 12 y 12 


12 


^ r 12 " y2 ) 


V2 


dy 


12 


-r 


12 


dx. 


0 0 A 12 J 12 ,_2 ,2 s 1 

?o - y?oV/2 ^ X 12 +y 12 


(fl2 ~ y 12j 


0. 


Similarly it can be proved that K = 0. 

X 31 y 31 


f -r 0 , r 2 r 


K. 


X 


23 y 23 


■23 


i r, 


23 J 


"31 


X 23 7 23 (2 r 12 + 2 r 31 )2, r o0 ^23 ^23 


-2 r 


12 


'31 • 23 


dx o:! dy. 


r 23 (r 12 ” r 31 ^ 


(2.V0) 


It is assumed that the shift of the roller is approximately 

same on both the sides from the mean position and hence 

K „ is approximately zero. Hence, the random varia- 
23 y 23 

bles x. . and y. . can be assumed to be uncorrelated. 

1 3 «3 



The mean and standard deviations of the various 


random variables can be derived as follows j 7 \ : 


The marginal density function, f v (y o0 ) , is given by 

•^23 2i 


f X 23 ^ 


f 2 r 31 

i ^ 

= I 2 lr 12 + P 31 > 2 r 2 - 3 ®23 > 

J" 2 r 12 


(2.42) 


f r 




and 


y 23 


2 

^r y 23 


23 


J-r 


23 


y 23 2 r 23 dy 23 


f r 23 


by. 


i -r. 


23 


'23 2 r 23 ^23 


0 , 


r 2 

- fai 


(2.43) 


(2.44) 


Considering the tolerances on the link lengths as normally 
distributed (with + 3<yband), the mean and standard devia- 
tions of the various random variables are given below : 


Random Variable 

1 • 

| A a 

Ab 

X 12 

*3-1. , 

y 23 

Mean 

\ 0 

0 

0 

0 

0 

Variance 

I (A a/3) 2 

J 

(Ab/3) 2 

x 2 2 A 


4/3 


From Eqs . (2.35) it can be observed that the 
acceleration and the Jerk are functions of the random 
variables A a, A b. x 12 , y 32 3X1(1 X 31 * Using the 

partial derivative rule, the standard deviations of the 
acceleration and the jerk can be computed as follows: 
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2 

) 








(2.13) 

( 2 . 14 ) 


where ( <TJ) . = standard deviation of acceleration at 

C* = C* 3 , 

(CT4) . = standard deviation of jerk at c* = o< . , 

J J 

A. = Acceleration at i= J . , J. = Jerk at ~\= V. 
J ^ ' U J d 

X = mean design vector = Null vector 

= ( 0 , 0 , 0 , 0 , 0 ) T . 


The partial derivatives, required for the compu- 
tation of standard deviation of acceleration and jerk, have 
been derived in closed form and are given in Appendix A. 


Optimum Allocation of Tolerances and Clearances : 

Very small values of tolerances and clearance, 
though give rise to small deviation in the acceleration 
and jerk, result in high manufacturing cost. Hence, maxi- 
mum possible values of tolerance and clearance should be 
allocated to various link lengths, keeping the deviation 
of the acceleration and jerk within tolerafcl® limits. This 
problem can be formulated as an optimization problem as 
follows : 
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Find X = 


which minimizes OBJ = 


5 


" x 1 N 


(A a. ' 

x 2 


1 ^ b 

1 X-. j 
( 

} J 

1 x i2 y 

x if ! 

| 

y 23 

x.- 

L j 


X 21 

t / 


X i 


subject to the following constraints 

1 i ^ x i ^ U i * i = 1 to 5 

c qpj ^ Ej 
c ap 3 ^ F d 

where 1^ and Ih denote the lower and the upper bounds 
on x, (i = 1 to 5) and E. and F. represent the maximum 

J. ti J 

permissible values of Q~£ and O’J at cX= OCy Tiie objec- 
tive function, OBJ, represents a measure of the manufac- 
turing cost. The constraint on the standard deviation of 
acceleration lias been incorporated to limit the magnitude of 
inertia forces. In a similar manner, to ensure smooth ope- 
ration of the mechanism, the rate of change of acceleration 
(that is jerk) has also been controlled by placing an 
upper bound on the standard deviation of the jerk. 


The nonlinear programming techniques (interior 
penalty function method coupled with variable metric 
method of unconstrained minimization and cubic Interpo** 
lation method of one dimensional minimization) have been 
used to solve this problem (6f . 
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NUMERICAL EXAMPLE: 

The data of the numerical example considered is 
given below. 

Upper bounds on the design variables : 

U ± = ( 3$ a, \% b, 0.2, 0.2, 0.2) 

Lower bounds on the design variables : 

1 ± = (0.2$ a, 0.2$ b, 0.01,0.01,0.01) 

Range of the input angle for a 4 slot mechanism 
with crank length of 10 cm is from - V?° to + 45°. The 
standard deviations of the acceleration and the jerk have 
been calculated at 5° interval. The values of E.. and 
are taken as : 

E. = 20$ Mean acceleration at (X= (Xjj 

3 

Fj = 30$ Mean j erk at (X= CX -j • 

The results of optimization are shown in Table 
2.7. There has been a reduction of 35*2$ in the objective 
function. The upper bounds on the standard deviation of 
jerk at o<= 0° andCX^ 5° have been found to be active at 
the opt imum point. The reduction in the objective function 
with the number of one dimensional minimization steps is 
shown in Fig. 2,10. The + 3 band widths of accelera- 
tion and jerk corresponding to the starting design vector 
are given in Table 2.f, . The band width corresponding to 
the optimum point are shown graphically in Figs. 2,11 and 
2.12. The band widths can be seen to be large at the 
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central part compared to those at the ends. The band 
widths are larger at the optimum point compared to those 
at the starting point. This reflects increase in toleran- 
ces and clearances at the optimum point compared to those 
at the initial point. 

2.5 PROBABILISTIC ERROR ANALYSIS CONSIDERING TOLERANCES 
AND CLEARANCES FOR THE INTERNAL 
GENEVA MECHANISM 

The probabilistic model of the internal Geneva 
mechanism is similar to that of the external Geneva mecha- 
nism. The Eqs. ( 2 . 13 ) and (2. 14 ), for the standard devia- 
tion of acceleration and jerk, apply here also. The partial 
derivatives, required for the computation of standard 
deviations, are given in Appendix A. The formulation of 
the optimization problem is the same as that for the exter- 
nal Geneva mechanism. 

NUMERICAL EXAMPLE: 

The data of the numerical example considered is 
given below. 

Upper bounds on the design variables : 

U ± = ( 3% a, b% b, 0.2, 0.2, 0.2) 

Lower bounds on the design variables : 
l jL = ( 0,2% a, 0,2% b, 0.01, 0.01, 0.01) 

Range of the input angle for a 4 slot mecha- 
nism with crank length of 10 cm. is from - 135° to 
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+ 135°. The standard, deviations of acceleration and jerk 
have been calculated at 5° interval and the upper bounds 
on the standard deviation of acceleration and jerk are taken 
as = 0.05 and F. = 0.4 for all j's. 

In this case it is not feasible to take E.. and 
Fj as fixed percentages of the mean acceleration and the 
mean jerk respectively (as has been done for external 
Geneva mechanism) because the mean values of acceleration 
and jerk are extremely small . 

The results of optimization are summarized in 
Table 2.9. It can be seen that the objective function has 
been reduced by 91 *3%. The upper bounds on the design 
variables a, Ab, x^ 2 > and have been found to 

be active at the optimum point. Among the behaviour cons- 
traints, the upper bound on the standard deviation of jerk 
at o(= 135° has been found to be critical at the final 
point. The progress of optimization with respect to the 
number of one dimensional steps is shown in Fig. 2.13* 

The ± 3 <3" band widths of acceleration and jerk are given 
in Table 2. 10 for the starting design vector and are shown 
graphically in Figs. 2.14 and 2.15 for the optimum design 
vector. It can be seen that the band widths are small at 
the middle compared to those at the starting and ending 
angles. The band widths can be seen to be much larger 
at liie optimum point compared to those at the starting 
point. This is to be expected since the starting design 
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vector is nearer to the lower bound vector while the 
optimum design vector is closer to the upper bound vector 


DISCUSSION OF THE BESULTS : 

The optimization programmes have been run on 
DEG 1090 system (a fourth generation computer). The opti- 
mization problem for the external Geneva mechanism took 
1.22 (1 min. and 22 sec.) CPU (Central Processing Unit) 
time while for the internal Geneva mechanism problem the 
time taken was I. 38 . From Tables 2. '7 and 2.9, it can be 
concluded that the objective function decreases rapidly 
in the first few iterations. If computer time is also a 
criteria then it is sufficient to optimize for two values 
of the penalty parameter i*. . For the subsequent values of 

K. 

more than two, the decrease in the value of the objec- 
tive function is not significant. Another interesting 
result is that for both the mechanisms, the active cons- 
traints at the optimum point correspond to upper bounds 
on some of the design variables and upper bounds on the 
standard deviation of jerk only. The standard deviations 
of the acceleration at various crank angles are well 
within the Units at the optimum point. 



31 



SUMMARY OF NUMERICAL RESULTS OF DETERMINISTIC ERROR ANALYSIS OF EXTERNAL GENEVA MECHANISM 
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TABLE 2-10 



















P = Pin on the driving member D 

D - Driving wheel (input wheel ) 2 f = Indexing period of the 

F= Follower (out put wheel ) driver 

L = Locking arm 


FIG. 2-1 .-a EXTERNAL GENEVA WHEEL AT ENGAGEMENT 




D = Driving wheel 

P = Pin (roller ) on the driving member 

F = Follower 
L = Locking arm 

2<f- Indexing period of the driver 

3F!G. 2 -l b INTERNAL GENEVA MECHANISM AT ENGAGEMENT 
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IG.2-1-C COMPARISON OF ACCELERATION .IN EXTERNAL AND 
INTERNAL GENEVA MECHANISM 
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at engagement 


Axis of slot at disengagement 
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CHAPTER 3 


CAM - FOLLOWER SYSTEM 

3.1 GENERAL DESCRIPTION: 

In this chapter, the follower responses of cam- 
follower systems is calculated probabilistically by consi- 
dering the errors in various design parameters for the 
following two cases ; 

1. Effect of errors in the geometrical parameters 
on the kinematic response. 

2. Effect of errors in geometrical and other system 
parameters (such as equivalent mass, stiffness 
and damping) on the dynamic response. 

For simplicity, the analysis is carried for a 
two dimensional disc cam with translating roller follower. 
The kinematic response has been calculated for rise period 
only while the dynamic response has been computed for both 
rise and dwell periods. 

The synthesis problems, involving the allocation 
of tolerances on the system parameters have also been 
solved using nonlinear programming techniques. 

3.2 KINEMATIC RESPONSE OF THE FOLLOWER: 

Figure 3*1 shows a disc cam with translating 
roller follower. The nomenclature is given below: 
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r 

r b 

a 


F 

F 

0 


1 

2 

2 


9 

^2 


S 1 


L 


o 


-f 


0 

*0 

s* 

1 


=* roller radius (cm) 

= radius of the lase circle of the cam (cm» ) 

= offset or eccentricity (cm.) 

- coordinate system fixed to the follower (X^ , ) 

= coordinate system fixed to the cam disc (Xg, ^ 2 ^ 
= cam rotation angle 

= parameter of the follower surface 
= constant angular velocity of the cam shaft 
( radians/ secs. ) 

= displacement of the follower (cm.) 

= total lift of the follower (cm.) 

T 

= design vector (x^ , x 2 , x^, x^, x^) 

corresponding to the tolerances in the geome- 

T 

trical parameters ( Ar, Ar^, /-a, AX, AY) 

= cam angle for total rise 
d s 1 

- ary 

o o 1/2 

= (r£ - a 2 ) (3.1) 


X 2 , ^2 = coordinates of the point of contact P. 

The condition of contact can be derived as 
follows ] 9 1 i 


E Sin 9 + F Cos 9 + G = 0 (3.2) 

where 

E = - a , 


and G 


0 
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From this equation, the parameter 6 can be rea dil y 
evaluated as 

i s i + 

© = tan" 1 ( -3 ° ) (3.3) 

sJj - a 

The coordinates of the cam profile can be derived 

as ) 9 1 J 

X 2 = - r Cos (© + 0 2 ) + a Cos 0 2 + (s^ + s c ) Sin 0 2 

(3.4) 

Y 2 = - r Sin (© + 0 2 ) - a Sin 0 2 + (s^ + s Q ) Cos 0 2 

(3.5) 

Using these equations, the derived relation for 
follower response for error analysis can be derived as 

s 1 = X 2 Sin 0 2 + Y 2 Cos 0 2 - r sin © - s Q (3.6) 

Equation- (3.6) indicates that s 1 , the follower 
response, is not an explicit function of the geometrical 
parameters, because the variables ©, X 2 and themselves 
are functions of s 1 as given by Eqs. (3.g) to (3.5), 
respectively. 

The tolerances in the geometrical lengths 1 r> and 
' r b ' are, obviously, random variables. Further some error 
is bound to occur in fixing the cam follower system and 
hence the offset or eccentricity ’a' can be assumed as a 
random variable. In manufacturing the cam, errors are 
unavoidable. Thus, the tolerances in the coordinates of 
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the cam profile can he assumed as random variables. Though, 
these coordinates change with the cam rotation angle, the 
tolerance ^ X in the X 2 coordinate and AX in the Y 2 
coordinate are supposed to be the same for all the cam 
rotation angles. It is assumed here that the workmanship 
of the manufacturer remains the same. 


To start with, a theoretical s^ is assumed and 
substituted in Eqs. (3.3) to (3*5) to get ©, X 2 and Y 2 
respectively. These values are then substituted in Eq, 
(3*6) to get s^ . The standard deviation of the kinematic 
response can be obtained by applying the partial 

derivative rule as 


( < r s 1 )j 


i- ( I 

1=1 V ^ X i 7 1 x x 3 


1/2 


(3.7) 


where 


( Va 

(s l)j 


standard deviation of s^ at 0 2 


s 1 at 0 2 = (0 2 ). 


^ 2^3 


X = mean design vector 

= (0.0, 0.0, 0.0, 0.0, 0.0) 


T 


The partial derivatives, 



are given below: 


required for the computation of 


?>s 1 

'^s 1 


Sln0 2 
Cos 0 2 
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■} s 1 

> x 3 

^ S 1 

3*4 


= - Sin 9 


y r n r; 9 3 S o 

- T Cos 9 . - -j-jr 


d S 1 

3 x 5 

3_9 

-3r b 

a q . 

3 a 
3 s o 


3 r b 


- tCose f!-irf 

v {<■; 

<S> ( v vo; 


„ / (y£ q 2 n 1/2 

r b / (r-b - a ) 


a) + (s-j + s Q ) 


: 1 


a) + (s 1 + s Q ) 


2 

2i 


7)S 0 

fa 


— / / 2 2 \ 1 /2 
= - a / (r b - a ) 


' PROCEDURE TO FIND THE MEM AND STANDARD DEVIATION 
OF THE KINEMATIC RESPONSE 

The various steps involved in the computation of 
the mean and the standard deviation of the follower response 
are given below: 

1. The dam angle for total lift, that is 0 q , is divided 
into a convenient finite number of intervals (n - 1). 

The discrete values of 0 ^ are numbered as (j^)^ * 

3 ’ = 1 , 2, .. . n , Set j=1, 

2. A suitable lift curve is assumed. Here cycloidal 
motion has been selected which gives motion with much 
reduced discontinuity of acceleration compared with 
other common types of motions. The lift s-j , 



corresponding to any cam angle 0 2 , according to cycloi- 
dal motion is given by 


L 

S 1 = T\ 

and hence 

, d s 

S 1 = 


f ~ ( I } Sin (27V 0 2 /0 o ) 


(3*9) 


1 


= ( 


L 

0 \ 

J £ 5 


r 


1 - Cos (2 7T 0 2 / 0 q ) 


1 


(3-10) 


The values of s. 


and 


particular value of 0 , 


= ( 0 2 ) 


tituted in Eq. (3*3) to get the 
9 at 0 2 = 0 2 _. (that is 9.. ) . 


are calculated for a 

. and are then subs- 
J 

corresponding value of 


3. The mean lift (s.). and standard deviation ( o~~ ). 

• 3 s.^3 

are calculated using Eqs. ( 3 . 6 ) and ( 3 * 7 ) respectively. 

4. Set new 3 = present 3 + 1 

and repeat steps 2 through 5 " until the 
angle ( 0 2 ) n is exhausted. 


3.3 OPTIMUM ALLOCATION OF TOLERANCES 

In most of the mechanisms that use cam follower 
systems, the follower response is not allowed to deviate 
more than a prescribed value from the mean (or nominal) 
position. The motto of the optimization problem is to 
allocate maximum tolerances to the various geometrical 
parameters (and hence minimize the manufacturing cost) 
keeping the follower response within the tolerable limits. 
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The tolerances for all the geometrical parameters 
are assumed to be normally distributed with + 3 band. 
Hence, if the mean value of one of the parameters is r 
and its standard deviation is ■ r - T , then the allowable value 
will be <r + 3 T- r ). Hence cr^, = Ar /3 = *3/3 where 
J\ r is the tolerance given to r. 


Similarly 

ffrb = *V3 = V 3 

(f~ a = Aa / 3 = x j /3 

'T"x 2 = AX/3 = x/3 

^2 = *V3 = x 2 /3 

The optimization problem is formulated as 

follows : 

Find T = ( X 1t x 2 , x 3 , x/ 

which minimizes 

5 

OBJ = ^ — L 

i =1 x i 


Subject to the following constraints 

X i ^ x i ^ u i » i = 1 to 5 

( cr~ ). 4 e 

s 1 j ^ 

where 

1 ^ = lower bound on x. , 

u^ = upper bound on x^ } and 
E = upper bound on the standard deviation of 
the follower response. 
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NUMERICAL EXAMPLE 

For solving the optimization problem the following 
data has been utilised. 



= (0.001, 

0.001, 0.3$ r, 

0.3$ r b , 

0.3$ 

a) 

{ u±, i 

•s 

CM 

• 

O 
' 

II 

0.2 , 3.0$ r, 

3.0$ r b , 

3.0$ 

a) 

r = 

1 .0, 

r b = 2.5, a 

= 1 .2 




Total angle of rise 0 Q - 150°. 


The problem has been solved for four different 
values of E, namely, .01, .06, 0.10 and 0.20. 

The results of optimization are summarized in 
Tables 3.1 to 3.*+. It can be seen that the objective func- 
tion has been reduced by 70 . 5 % for E = 0.01, 90 . 7 % for 
E = 0.06, 91.5$ for E = 0.1 and 91.8$ for E = 0.2. The 
progress of optimization for E = 0.1 is shown in Fig. 3.2. 
The objective function decreases rapidly in the first few 
iterations for all E’s except for E = 0.01. The standard 
deviation of the response corresponding to various values 
of $2 "^* ie optimum point are given in Table 3.5 for 

different values of E. It is important to note that the 
standard deviation almost remains the same for all values 
of particular X. This is the reason why at 

the optimum point corresponding to E = 0.01, all the 
upper bounds on the standard deviation of lift are critical. 



For E = .06, E = 0.1 and E = 0.2, none of the constraints 
corresponding to the upper bounds on the standard deviation 
of lift is critical, revealing that the upper bound has been 
relaxed too much. There is no significant improvement in the 
objective function at the optimum point by relaxing the upper 
bound on the standard deviation of lift from E = .06 to 
E = 0.2. In most of the practical cam follower systems, 
the value of E is to be restricted to be less than or equal 
to 0.01 . 

3.4 DYNAMIC RESPONSE OF THE CAM FOLLOWER SYSTEM 

The kinematic analysis alone will not be sufficient 
for high speed cam mechanisms. It needs to be supplemented 
by the dynamic analysis to ensure the faithful reproduction 
of the follower motion. As the exact mathematical treatment 
is extremely complicated, mathematical modelling of the 
whole mechanism is necessary. Assuming the cam shaft to be 
absolutely rigid, the system can be idealized as a two degree 
of freedom system lisp By reducing the load and the entire 
follower train into a single mass m, the system can be 
further simplified so that it has a single degree of freedom 
as shown in Figure 3.3.b. 

The objective of the dynamic analysis of the cam 
system is to determine the response x when s-j is known 
as a function of time. 



The nomenclature used is given below: 
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ii 


C 1 

c 2 

N 

C*~ 


o 

S 1 


- Equivalent mass of the follower train (Kg) 

= Stiffness of the follower spring (Kg/ cm. ) 

= Stiffness of the follower train (Kg/cm.) 

= Damping coefficient of follower (Kg-sec,/cm. ) 

~ damping coefficient of follower train (Kg-sec./cm) 
= Speed of the cam shaft (rpm) 

= Angular velocity of cam shaft = — radians/sec. 

= Kinematic response of the cam (cm) 
d s 1 

= “dF" 


x 

— *- 
XX 


(<ri ) 

x i " 





- Dynamic response of the cam follower system (cm.) 

= (r, r b , a, m, k 1 , k 2 , c 1 , c 2 , X 2 , Y^) 

= Standard deviation of lift during the rise 
period at the ith time step 

= Standard deviation of lift during the dwell 
period at the ith time step 
* Design vector = ( £ r, Aa, Am, £ ^ , 

Ak 2 , Ac^ Ac 2 , Ax, £Y) 

- Vector of upper bounds on the variables X 
= Vector of lower bounds on the variables X* 

= Upper bound on the standard deviation of 

lift (cm) 



Numerical Solution for Dynamic Response 
During the Rise Period ~ 
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The dynamic equilibrium equation of the system 
shown in Fig. 3.3»h is given by 


m x + (c-, + c 0 ) x + (k- + k 0 ) x 


l^s^ + C2 S-j 


or 


c. + c 0 . k., + k 0 

x + ( -J — — 2 ) x + ( -3 L- — 2 ) x = ( 


k2 Sl *f" ^2 


m 


m 


m 


) 


(3*11) 


The right hand side of the above second order 
linear differential equation, with constant coefficients, 
is time dependent. Here the kinematic response (s^) is 
given by the following expression 


S 1 = X 2 Sin ^2 + ^2 Cos ^2 " r Sin G “ s 0 (3.12) 


It is not possible to get the time derivative 
of s^ explicitly (as the time derivatives of X 2 and Y^ 
are not known) and hence the dynamic response x can not 
be obtained as a closed form expression. In this work 
*s^ is calculated by numerical differentiation. 


Equation (3.11) can be written as 


x+2nx+p^x = q^ 


(3.13) 


where 


n = 


C 1 + c 2 
2m 


2 k 1 + k 2 


m 


q i = 




m 


and 


forcing function per unit mass 
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In most practical cam mechanisms the values of n and 
p are such that the roots of the auxiliary equation 

• • *2 

x + 2 n x + p = 0 are complex and hence the 

system under consideration is underdamped. The transient 
response is given by 


-nt 


x 


= e (A.j Cos p d t + A^ Sin p d t) 


( 3 . 1 ^) 


where 


f 2 ..a - ! 1/2 


p d = j p - n 

To determine the constants A^ and A^, it is 
assumed that at the initial instant (t = 0), the vibra- 
ting body is displaced from its position of equilibrium 
by the amount x Q and has an initial velocity x Q , 


Hence 


A. = x^ and A 0 
1 o 2 


x + n x 
o o 


p d 


Hence transient response = e~ n ^ (x q cos p d t + 


x + n x„ 
o o 


p d 


Sin p d t ) 


Tho response for an arbitrary disturbing force 
q = f(t’)is given by the Duhamel's integral as 


-nt r 


X = 


P, 


) 


nt 1 


q . Sin p d (t - t* ) dt f 


The complete solution is given bye the following 
expression : 
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x = 


4* 

Pd 


r° 

.t 


o 


■‘■CL 
nt 1 

e q Sin p d 


u. 


1 

(t - t’ ) dt*J 


(3.15) 


As a particular case when the forcing function is a step 
function, as shown in Fig. 3.^+* the response is given by 


x = 


t- 

*1 


nt' ”1 

+ ^ J Q e Sin Pd (t ~f) dtj 


Ql 

where q = — - . 

H m 

The integral in the above equation is integrated by parts 
to get 


x 




Cos p d 


*fc 4* 



+ n x o 

Pd 



4 




r nt (Cos p d t + ~ sin p d t) 


I 


The usual procedure is to approximate the forcing 
function by a series of step functions as shown in Fig. 3. 4. 
Care is taken to see that the time step is sufficiently 
small compared to the natural time period of the system. 

In axqr time interval t d _^ ^.t-^t^, the response of the 
one degree system may be calculated as the sum of the 
effects of the initial conditions at time t^ ^ and the 
effect of the impulse within the interval /% t^, as 
f ollows : 
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-n (t - t. ,) p 

X a: e i x i _ 1 Cos P d (t - t^) 


x . ' j + n X. . 

+ Q± n (t - 


Pd 


q. ( -n (t - t. 1 ) 

+ -4 < 1 - e :L ~ 1 

P 


Cos p d (t - \_'j) 


+ ^ Sin Pd (t - *1-1 } J j 

At the end of the interval this expression becomes 


x^ = e 


-n A t^ 


x i-1 Cos p d & t ± + 

P -n „rixt.. / 
it- e 


x ±-i + n *i-i 


Sta p d AtJ 


1 (Cos P 4 A *i P an P d A ^ 


(3.16) 


The velocity can be determined from the following expre' 
ssion 





Pd 


-n At d 


x i-1 + n 'S.-l 

- *1-1 Sin Pa Ah + ' * Pd 1 


n X i-1 + n X i-1 

- (x i-i Cos Pd A % + ~ Sin Pd 


q. - n A t . „2 

+ * e ( 1 + \ ) Sin p d /t. 

P P. 


Cos P d At ± 


4*1 > 
(3.17) 


Equations (3.16) and (3*17) represent the 
recurrence formulas for calculating the dynamic response 
at the end of the ith time step. They also provide the 
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initial conditions of displacement and velocity at the 

"t In 

beginning of (i + 1) step. Continuing this way, the 
displacement x (0 Q ) and velocity x (0 ) are found. 

Analytical Solution for Dynamic Besoonse 
During the Dwell Period 

During the dwell period the rise of the cam is 
constant and is equal to L . Substituting s^ = L Q and 
3^ = 0 in Eq. (3.11) the dynamic equilibrium equation 

modifies to 


• • c. + c 0 • 

x + ( x + ( 


k. + k 0 k~ L 

1 2 ) X = 2 


m 


m 


(3.18) 


The solution of this second order linear differential 
equation with constant coefficient gives 


x = e (A 3 Cos p d t + A^ Sin p d t) + ^ + k ' 2 

The initial conditions used to calculate the constants 
A^ and , are as f ollows : 

0 yj 

At t = , x = x (0 ) 

180 o.' ° 

0 7T 

At t = — — - x = x (0 ) 

180 CO ° 

The constants can be evaluated as 

k p L 

A 3 = x - k 1 + t 2 

\ = <* (0 O > + n A 3 ) / p a 
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Calculation of Standard Deviation of Response 


The partial derivatives required for the calcula- 
tion of the standard deviation of lift during the rise 
period are calculated by numerical differentiation. 


The standard deviation is given by 

f 10 -x. 

(<J x i ) r = ( £xx(k) } 


c x i 


XX(k) f 


( 3 . 19 ) 


The partial derivatives needed for the calcula- 
tion of standard deviation of lift during the dwell 
period are calculated analytically. It can be seen from 
equation (3.18) that the dynamic response in the dwell 
period is independent of the geometrical parameters. 






7b 


Calculation of Standard Deviation of Respon se 


The partial derivatives required for the calcula- 
tion of the standard deviation of lift during the rise 
period are calcinated by numerical differentiation. 

The standard deviation is given by 


r 




10 

y" 

s 

k=1 


"8 X i 


S'XX(k) 


2___ v 1 

XX(k) J 


1/2 


(3.19) 


The partial derivatives needed for the calcula- 
tion of standard deviation of lift during the dwell 
period are calculated analytically. It can be seen from 
equation (3.18) that the dynamic response in the dwell 
period is independent of the geometrical parameters. 
Hence 


(<r~ ) 

x i d 


8 

£ 


L k=4 


3 x i 

d'XX(k) 


2 

) 


2 _^ 

0^ XX(k)| 


1/2 


( 3 . 20 ) 


where 

^ x i 
^XX(k') 


-nt _ p d 

e (- A 3 Sin p d t . p d £xx(k) 

+ ^ COS p a t . pd 


-nt n 

- e n - T&Ck ' J (a 3 Cos *' :+ \ Sin p a 


CON = 


k 2 L o 
k 1 + k 2 » 


[ on n 3 n 

3 XX (5) - ~ m ’ IWI 

8 n fen 1 

535rt7T = JWm ~ 2m 


^n 

7) XX(6) 


= 0 
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3.5 ALLOCATION OF OPTIMUM TOLERANCES: 

The tolerances in the geometrical parameters and 
the variations in the system parameters are assumed l 

to he normally distributed. The optimization problem is 
formulated as follows: 

Find X = (-'Nr, Ar c , *\a, Am, .^k 1 , ^>k 2 , r\c 1 , 

zic 2 , -^X, ^-Y) 

which minimizes 

1 ° ’ . 

= X rizr 

ill 

Subject to the following constraints 

L (I) <. X (I) ^ U (I) , 1 = 1, 10 

( (SV) E, i = 1, 2, ... n. 

where L (I) = Lower bound on X (I) random variable 
U (I) = Upper bound on X (I) random variable 
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( 6"x^) = Standard, deviation at ith step 

E = Upper bound on the standard deviation of 
the dynamic response x. 

NUMERICAL EXAMPLE 

For solving the optimization problem the following 
data has been utilised. 


U(I) 

= 3 % XX(I) , 

I = 

1, 3 

U(I) 

= 10$ XX(I) , 

I = 

OD 

J- 

U(9) 

= U(10) = 0.2 



L(I) 

= 0.2$ XX(I) , 

I = 

1, 3 

LCD 

= 2 % XX (I) , 

I = 

00 

L(9) 

= L(10) = 0.001 




Rise period = 

0 to 

i5o° 


Dwell period = 

150° 

to 210 


XX(I) = ( r, r b , a, m, , k 2 , c^, c 2 ) 

= (1.0, 2.5, 1.2, 0.3, 40.0, 700.0, 0.005, 

0.005) 

N = 3000 r.p.m. 

Figure 3.5 shows the mean dynamic response, the assumed 
motion and the follower command for various values of 
0 2 . 

The optimization problem has been solved for 
three different values of E, namely, 0.06, 0.1 and 0.2. 
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The results of optimization are summarized in 
Tables 3.6 to 3.8. It can be seen that the objective 
function has been reduced by 26. for E = 0.06, 57 . 0 % 
for E = 0.1 and 66 . 8 % for E = 0.2. Once again the objec- 
tive function decreases rapidly in the first few iterations. 
The progress of optimization for all E* s is shown in Fig. 
3.6. The standard deviations corresponding to the various 
values of 0 ^ at the optimum point for different values of 
E are given in Table 3 . 9 . It is to be noted that the 
standard deviation increases with the increase in the 
dynamic response up to the rise period. From Tables 3.6 to 
3.8 it can be observed that as the value; of E decreases 
from 0.2 to . 06 , more and more constraints corresponding 
to the upper bound on the standard deviation of lift 
become critical. From Table 3«9 ? It is possible to plot 
± 3 6 ~ band about the mean dynamic response for any value 
of E. 

It can be concluded by observing the standard 
deviation of the kinematic response given in Table 3.5 
and the standard deviation of the dynamic response given 
in Table 3»9j that the former is much smaller compared 
to the later. Thus the consideration of kinematic analysis 
alone is not sufficient to synthesize the system. It has 
to be supplemented by the dynamic analysis which gives 
more realistic results. 
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TABLE 3.8 

Initial point x-, = (.003, .006, .003, .00001,1.0,20.0, .0002, .0002, .005, .005) 
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PF = Penalty function 
0 B J = Objective function 

ITER = Cumulative number of one dimensional 
minimization iterations 
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FIG. 3-2 ITER 

PROGESS OF OPTIMIZATION FOR CAM PROBLEM 
(Kinematic response ) 
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PF= Penalty function 
OBJ = Objective function 

E = Upper bound on the standard deviation of 
dynamic response 

ITER = Cumulative number of one dimensional 
minimization iterations 
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FIG. 3-6 ITER 

PROGRESS OF OPTIMIZATION FOR CAM PROBLEN 
(Dynamic response ) 



CHAPTER ^ 


GEAR TRAIN 

4.1 PRECISION GEARING: 

Gear applications vary so much that the design 
requirements range from minimum controls and maximum 
loose tolerances, to complete controls and minimum tole- 
rances of the order of 0.0001 inch. Gears of high quality 
with tolerances in the order of thousandths and ten- 
thousandths of an inch, are called precision gears. These 
occur in all types, pitches, sizes and materials. 

Precision gearing in its applications is largely limited 
to instrumentation, automatic controls (including servo- 
mechanisms), computers, machine tools and special high- 
speed/or high-power units. 

A relatively recent innovation is the specifica- 
tion of stringent performance requirements upon gearing as 
one component in a complex system. Many precision gears 
used in the electronics industry and military weapon 
systems are rigorously assessed for performance and error 
contribution. The art of gear technology has been trans- 
formed to an engineering science encompassing rigorous 
mathematical analysis and error predicting techniques. 

The primary objective of a precision gear is to 
provide the functions of motion, power transmission and 
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dynamic response with a high degree of accuracy. In order 
to achieve the desired level of performance, the various 
sources of error within the gear train have to be controlled. 
The two major errors in a gear train are transmission error 
and backlash. The overall precision of a gear train is 
often specified by a quantity known as integrated gear 
train position error, which is the net sum of transmission 
error and backlash. These errors and their sources are 
defined in Appendix B. The probabilistic analysis of the 
various errors was given by Michalec j 5 j and is summarized 
in Appendix B. With the availability of nonliner optimiza- 
tion methods, it is possible to optimally allocate the 
various tolerances (errors) in a gear train while satis- 
fying the specified requirements of precision. 


if .2 ALLOCATION OF ERRORS IN A GEAR TRAIN 

In general synthesis problem can be stated as 
an optimization problem as follows : 

^ T 

Find X = (x 1 x 2 x^ x^ , . . . x n ) 

n 

which minimizes OBJ* = / _ — 

i=1 i 

and satisfies the constraints 
1 i ^ x i ^ u ± 

/^I ^ ^I^max. 

'^""i^max. 
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where 


x_p i = 1, 2, . .., x = errors (tolerances) to 

he allocated, 

1^, u^ = lower and upper hounds on x^, 

and (T^ = mean and standard deviations of 

integrated gear train position error, 
and 


^I^max. 311(1 ( 


maximum permissihle values 
of U i and CT^-. 


NUMEEICAL EXAMPLE : 


The gear train shorn in Fig. 4-.1 is considered for 


the optimal allocation of errors. The design vector is 
taken as 


t 

'x 1 
X 1 


maximum position error for gears of mesh 1 


x 2 


"ball hearing eccentricity for mesh 1 


x 3 


; maximum position error for gears of mesh 2 

i ! 

X = 

i 


= / 

> ^ 

! i 

maximum total composite error for gears L 

of mexh 3 

! 

x 5 

centre distance tolerance for mesh 3 


x 6 


gear size tolerance for gears of mesh 3 


*7 

^ 4* 


< component shaft runout for mesh 3 J 


The starting value, the lower and the upper hounds on the 
design variables are taken as follows (inches) : 
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i 

T 

i 


— r 

u i 

| (x^) starting 

1 


0.00005 


0.00100 

0.00010 

2 


0.00005 


0.00100 

0.00015 

3 


0.00005 


0.00100 

0.00015 



0.00010 


0.00100 

0.00050 

5 


0.00010 


0.00100 

0.00050 

6 


0.00050 


0.00200 

0.00100 

7 


0.00010 


0.00100 

0.00050 


It is assumed that the integrated gear train 
position error between the encoder shaft S-h and the out- 
put shaft S-1 is to be limited as 

( = 30 seconds of an arc and 

( CTI) v = 10 seconds of an arc. 

JL TClQjjL 

The other pertinent data of the problem is given 
in Table .1 . 

COMPUTATION OF jJt ^ and CPJ : 

The mean and standard deviations of the integra- 
ted gear train position error ( j and CPj) can be 
computed by using the formulas given in Appendix B as 
follows : 

(a) Transmission error of mesh 1 s 

Error Sources : 1 Maximum position error 


2 Ball bearing eccentricity 
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The mean and standard deviations of individual gear and 
pinion transmission error'" are calculated using Eq. (B.1) 
and those of mesh 1 using Eq. (B.2). The mesh transmission 
error is converted into angular measures as + 


ii 1 

/Vm-1 =2 ( >*T ) m -1 x 


1 / 180 x 60 x 60 


) secs 


( ^ I^m-1 


- g < ^T^m-1 s 4 


x i- ( 12Q- X-6Q 2L.60 ) secs. 

R-i 7V 


where = pitch radius of the gear in mesh 1 

= 2 . 1 87 5 inch . 


(b) Backlash of mesh 1 


The backlash error will not be there as spring- 

loaded antibacklash gears are used in mesh 1. Hence 

(a) (a) 

^B^m-l = ( ^"Vm-I = 0 

(c) Transmission error of mesh 2t 

Error sources : 1 . Maximum pisition error 

2. Ball bearing eccentricity 

3. Clearance between gear bore and 
shaft 

b. Shaft runout at gear mounting 




99 


Equations (B.1) and (B.2) are used to calculate the mean 
( T ) m _2 and standard deviation ( mesil 2 * 

These are converted to angular measures as follows : 


, (a) 

( J^T^m-2 

(a) 

^ '■^"~T^m~2 


where R r 


= 5 ( /lV2 x fc x ( 


180 x 60 X 60 

7T 


secs, 


_ 1 { ^ -.^1 ^ / 180 X 60 X 60 \ 

- 2 ( ^ 1V2 x x ( > secs ' 

= pitch radius of the gear = 1.250 inches. 


(d) Backlash of mesh 2: 

As the gears are spring-loaded, 
(a) (a) 

( B^m-2 = ( = 0 

(e) Transmission error of mesh 3* 


Error sources: 1. 

2 . 

3. 


Total composite error 
Clearance between gear bore and 
shaft 

Shaft runout at gear mounting 
Ball bearing eccentricity 


After determining C and ( CTtj-,) m _ 3 , th e f °l lowin & 

formulas are used to find their angular measures : 


180 x 60 x 60 

’ TV 


) secs. 


1 1 1 - , (a) 

(^T)m-3 = 2 x ^/*aPm-3 X S3 ( 

( rr- - 1 T f n— _ 1_ ( 180 X 60 X 60 \ 

^ ^T^m-3 ~ 2 ( cr T ) m-3 H 3 ' TV 

where the pitch circle radius (R3 ) is 1.250 inch. 
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(f) Backlash of mesh 3“ 

Error sources: 1. Centre distance tolerance 

2. Gear size tolerance 

3. Gear size allowance 

4. Ball bearing eccentricity 

7. Clearance between gear bore and 
shaft 

6. Clearance between bearing outside 
diameter and housing bore 

7. Clearance between shaft and bear- 
ing bore 

The mean and standard deviations of backlash of mesh 3 are 
calculated according to Eqs. (B3+) and (B.5) . These values 
are converted to angular measures as 


(a) 

( M B^m-3 


1 

(180 

* 60 ^ 60 ) secs. 

E 3 

(a) 

( ^m-3 

= ( °Ts>bh3 X 

1 

(180 

x 60 x 60 N 
— ) secs. 

/ \ 

E 3 

where 

= pitch radius 

of 

gear 

= 1 .25 inch. 


(g) Integrated gear train position error: 

The transmission error and backlash between 
shafts S-1 and S-4 are calculated using Eqs. (B.3) and 
(B.6) with the velocity ratios to shaft 1 (V^ = 1, 

= 10, Y^ = 100). Finals, the mean and standard 
deviations of the integrated gear train position error 
are computed with the help of Eq. (B.7). 
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NUMERICAL RESULTS 

The results of optimization are summarised in 
Table 4,1. It can be seen that the objective function has 
been reduced by 32 . 6 $. The upper bounds on the design 
variables, (centre distance for mesh 3) and 
(component shaft runout for mesh 3) have been found to be 
active at the optimum point. Among the behaviour cons- 
traints, the upper bound on the mean value of the integra- 
ted gear train position error has been found to be critical 
at the final point. The progress of optimization with 
respect to the number of one dimensional steps is shown 
in Fig . 4.2. 
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TABLE 4.1 

Design Data for the Gea_ Train Shown ir Eig. 4.1 



Item 

7 

0 


Mesh Number 





r : 3 " 

5 1 5 

2 

m 

H Ir 

JL_? 

1 . 

Centre distance 
tolerance 

- 

- 

+ 

- 0 



2. 

Gear size tolerance 
(testing radius) 

- 

+ 0 
- x 6 

+ 0 

-0.0005 

+ 0 

-0.0005 

3. 

Gear size allowance - 
(testing radius) 

- 

0.0002 

0.0002 

0.0002 

4. 

Ball bearing (inner 
-race) eccentricity x^ 

x 2 


0.0002 

0.0002 


Outer race eccentri- 
city 

— 

0.0002 

0.0002 

0.0002 


Radial play 

*- 

•* 

0.0003 to 0.0003 to 0.0003 ' 
0.0006 0.0006 0.0006 
(>t= 0.00022. 

<3T= 0.00005) 

5. 

Maximum position 
error tolerance 

± X 1 

± x 3 

— 

- 

- 

6. 

Maximum total l 
composite error 

- 

- 

x 4 

0.0005 

0.0005 

7. 

Shaft run out at 
gear mounting : 

gear 

pinion 

— — * 

0.0001 

0.0001 

See item 
11 

See item 
11 

0.0001 

0.0001 

See item 
11 


Table contd. on next 
page 
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TABLE 4- .1 (Continued) 


Item | 

Mesh Number 



r 

1 5 2 0 3 

0 j, 

5 * 

5 tr 

0 5 

8. Clearance between 

Selected 



gear bore and 

for 0.0001 



shaft : 

to 0.0002 



bore diameter 
tolerance 

Clearance , n nnn 0 

/*= o.oooi? +°- 0002 
<= = 0.00002 “ u 

0.0003 

0.0003 

shaft diameter 

- 

0.0002 

0.0002 

tolerance 

Gear 
+ 0 

-0.0001 
Pinion 
+ 0 

-0.0002 



allowance 

0.0001 

0.0001 

0.0001 


9. Clearance between 
shaft and bearing 
bore 


0.0001 to O.OOOltD 0.0001 to 

0.0002 0.0002 0.0002 

(/ 1=0.000075 

T-=0. 00001 ) 


10. Clearance between 
outside diameter 
and housing bore 


0.0001 to 0.0001 to 0.0001 to 

0.0002 0.0002 0.0002 

( ;J.=0, 000075 
0 - = 0 . 00001 ) 


11. Compunent's tole- 
rance ( S-4 and 
S-6) 

(i) Shaft: 

Shaft run out 

Shaft radial play 

Shaft diameter 
tolerance 

(ii) Mounting: 

Allowance 


x* maximum 
0.0003 to 0.0006 

0 to - 0.0002 

0.0001 


20 ° 


12. Pressure angle for 
all gears (0) : 



Active constraints are as follows 

GG(5) = -0.07465207 = Upper bound on 
G6(7) = -0.09496765 = Upper bound on 

GG (15)=—0. 0001 391 545 =Upper bound on the mean of the total error 
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CHAPTER 5 


CONCLUSIONS AND RECOMMENDATIONS 
5.1 CONCLUSIONS; 

By considering the analysis and synthesis of 
mechanical error in Geneva mechanisms, cam-follower system 
and gear train, almost all types of problems in the field 
of mechanisms have been covered. The probabilistic approach 
used in this thesis can be considered as a more realistic 
procedure since, in practice, most of the design parameters 
are really random in nature. Although normal distribution 
has been assumed in this work wherever necessary, the basic 
approach will not change even if the variables follow some 
other type of distribution. 

The tolerance? in the centre distance (A a) 
and in the crank length (A 13 ) mainly govern the deviation 
of acceleration and jerk from the mean values in the external 
Geneva mechanism while the clearance in the slot width ( c) 
controls the change in acceleration and jerk about the mean 
values in the internal Geneva mechanism. The dimensions 
(a and b) and the number of slots are same in both the 
Geneva mechanisms. The initial feasible point is also same. 
At the optimum point, the objective function gets reduced 
by 35 . 2 % in the external Geneva mechanism and by 91*3$ in 
the internal Geneva mechanism. So if one of the two 
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mechanisms is to be used it is cheaper to use the internal 
Geneva mechanism. In this case and also in other optimiza- 
tion problems solved in this thesis, the major reduction of 
the objective function occurs in the first two values of the 
penalty parameter. This shows that if the computer time 
available is limited, it is advisable to run the optimization 
programme for only first two values of the penalty parameter 
and to take the point obtained as an approximation to the 
true optimum point. 

None of the authors in the past combined the 
kinematic and the dynamic behaviour in the error analysis 
and synthesis of cam-follower systems. In this work, the 
response of the follower has been calculated by considering 
the tolerances in the cam profile, tolerances in the geome- 
trical parameters and errors in other system parameters, 
like equivalent mass, stiffness and damping during the rise 
period. It has been revealed that there is no significant 
improvement in the objective function at the optimum point 
by relaxing the upper bound on the standard deviation of 
the kinematic response from 0.06 to 0.2. On the contrary 
the objective function improves considerably, namely from 
26.1$ (reduction) to 66.8$ (reduction), as the upper bound 
on the standard deviation of the dynamic response is 
changed from 0.06 to 0.2. The optimum point, in the case 
of the kinematic response, is controlled by the upper 
bounds on the geometrical parameters while in the case of 
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of dynamic response, the equivalent mass mainly governs 
the optimum point. An important conclusion drawn is that 
the kinematic analysis alone is insufficient to define the 
follower response and that it has to be supplemented by 
the dynamic response in order to predict the lift of the 
follower with a greater degree of accuracy. 

In the case of the synthesis of gear train, the 
mean error is not zero as in the cases of Geneva mechanisms 
and cam follower systems. Here the optimum point is con- 
trolled by the upper bound on the mean value of the inte- 
grated gear train position error. 

9 . 2 BECOMMEHDATIONS : 

A more realistic error analysis of the Geneva 
mechanism problems can be made by considering the frictional 
force between the slot and the roller. 

While calculating the response during the dwell 
period of the cam follower system, the effect of errors in 
geometrical parameters was not considered. Hence the 
present work can be extended by considering the errors 
in geometrical parameters also, in the dwell period. 

It is recommended that the strength based 
error analysis be coupled with the present error analysis 
in the synthesis problem of gear-trains. 
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All the synthesis problems can be solved by 
posing them in an alternative format as follows : 




Find X which minimizes 

OBJ 


N 

— 


j=1 


(cr 0 ) d 


and satisfies the constraints 
n 


£ 

i=1 


1 


x i 


<C M 


and 1.^ u ± 

where (U~0). is the standard deviation of the output 
variable 0 at station j (j = 1 , 2, N) and M is an 

upper bound on the measure of manufacturing cost. 
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APPENDIX A 


PARTIAL DERIVATIVES OF ACCELERATION AND EERK OF 
GENEVA MECHANISM WITH RESPECT TO 3HE RANDOM 

VARIABLES 


Let the random variables be V (i = i to 5)* 


Then 

v 1 = X 

(1) = Aa, 

V 2 

= X (2) =2Ab 


v 3 =X 

(3) = x 12 , 

\ 

= X (W = y 23 


= X 

(5) = x 31 




The partial derivatives of acceleration and jerk are required 
to be calculated at the mean design vector. From Eq. (2*33) 
it can be seen that 

A = 

1 = 

The partial derivatives of the equivalent lengths with 
respect to the various design variables are given below 


u e 

^ a e 

- *v 3 - 


3 a 
w e 

?,V 2 

_ a e 

- b\ 

c 

- ~ 

a 

e 

A 

= 0 

C b e 

rq 

2 T e 

1 , 

a A 

5 V 1 

3 b e 

- inr 3 

,3 

b e 

A 

= 0 

■J c e 

Tv^ 




A c 

% 

c 


= 1 

J 

e 

-Jq 

c e 

~ :-v 2 

v ^ 

~ T 

e 

¥ 5 

= 0 


A (a e’ A- °e’" < > 


A (a + ^ a + X -|2 J ^ + - x ^ + s X23 > ^ ) 

J < a e > b e > c e j ^ ) 

J (a + Ei a + x 12 , b + -A b + x... , y p - , < ) 


(A.1 ) 



1 12 


The partial derivative of acceleration with respect to the 
random variable V^ is given below; 


A 

i V. 

1 


A, 


= < V > <7 


c v. T ; v. 

i a 


B i 1 -'' A i ^ B i 

) + ( — ) ( V rr + 1 A rr ) ! 


Y- 


dV, 


a V, 


■ f(A^ + ) (A^ + B.j ) + (A^ + B^ ) | 


o t i c* A^j 

\ f A * Tb N /A t *n A / 1 


OB 


+ ) 
i - a 
B, 


+ ^ ^ A i + B i ^ ^ A i + B i ^ ^ B' Vj 

+ ( ^ (A i + B i^ (A i + B -\' ) +c ,^v7 ) 

2 -i 1/2 


where V = |l - (A^ + B^ ) 


(A. 2) 


It can be observed that 


C 5Y~ ), 

c V 1 i X 


/ a a x 

( *V 


x 


and ( ) 

' V 2 


X 


( ) 
} 


X 




where X = mean design vector 

= null vector = (0 , 0, 0, 0, 0) 


(A.3) 


The partial derivative of jerk with respect to the random 
variable can be derived as follows; 


c! J 


v tr* m 

a A, r B. 


IV, 


( f ) ( Jv7 + 3v7> + .b^ca^b;’) 


<j* i> r q o t t 

( TST + ~rq } + ( ^ KA 1 + V CA 1 + B 1 ); 


a 

,1s B 


*- 


;w. 


ciB i 

d\ 


) + ( -4) 


v~ 


(A 1 + B 1 )(A 1 + B 1 ). 
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It x H 

£ i 'B 1 


3' A-j A B„ 


( . + . ) + Ui + B 1> ( . + - ) (Ai + Bi ) 

■X* MM — ■* 


, (3 A., , t it n.l , "a . i t. 

+ ( — v7 + -^Y7 )(A 1 + Ej) (A 1 + B-j )^j + ( v 3> (A 1 + B 1 ) 


« 2 


£ A 


C. ^ d B i . , ,3 3 A. 3 B. 

( rvT + "7v7 ) + ( -^5 )(A 1 + B 1 ) (A 1 + B 1 ) ( ^Y7 + XT ' 7 > 


+ ( (A n + B^ 2 (A* + B^) 3 (. 4-^1 + ) 

a T i t 3 ^ A 1 -7 B.. 

+ ( J } f 2 (A i +B i } (A i + V c inr + w > 


A. 


B< 


t 2 


3 A, 


D B, 


{ f C. CL. K. C- ^ 

+ 3 (A.j + ) (A^ + B^ ) ( y- + 3 y~ ) 


where 


V 


2 .,1/2 

1 - (A 1 + B^) | 


(A.lf) 

(A.5) 


As in the case of acceleration, the following relations hold 

good in the case of jerk also: 

, ! 

a j >L = (^f)L 

1 X ' v 3 X 

(A.6) 


( 


TT 


and C Srr~ ) 


bV, 


3 . X 

£$■ J > 


X = < *? 


X 


A. 1 EXTERNAL GENEVA MECHANISM: 


From Eq, (2.6), it can he seen that 


A„ 


b sin S 

6 . 


(a e + b e 


2 a h cos ./• ) 
e e 


172 


Hence (A^) 


b sin o£ 


(a 2 + b 2 - 2 ab cos. a. ) 


172 


X 


(A.7) 



tt 

B 1 

i * m 
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Thus, to calculate B 1 , B 1 , B 1 , A 1 , A 1 , A 1 and A 1 
at the mean design vector X, it is necessary to change 
a e to a, b g to h and c 0 to zero, in Eqs. (2.6) to 
Eq. (2.16). The partial derivatives of A^ , A*, A^', A>"'» , b! ( , 

« Hf 

and B^ with respect to the random variable TT are given 
below : 


d A. 


?! 


where 

Y 


i -3 b 

( 1 ) ( i — i 
K Y } K ^V. 


Sin ) - ( — U- ) PQ . b Sin 

O 0 


< a e + / 


2 a b Co s d ) 
0 0 


2 Y 

1/2 


(A. 8) 
(A. 9) 
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- f ~ 'T d p 
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e c V, 


e )7 i 
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- 2 Cos l (a + b 


1 a 


)7 


e 


(A. 10) 


c' B. 


V 


= 


I V Tf 


y- L? v i 

> a 


( a e “ B 0 Cos x ) 


r 0 Q 

e 


+ c ( 2 2 _ 

e k dY. f >V 


Cos c 


*> 1 


“ ( ~ TjT ^ c e ^ a e “ b e Cosc< ) • P( 3 


(A. 11) 
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)b 


e e dV ± 


— )+( np) (- f) . PQ 18 a 2 b 8 Sin 2 *. Cos* 


e e 


a 


+ ( -C 1 ) * 18 Sin 2 / Cos o< (2 a rp » b 8 
Y ? e c v ± e 

+ 3 tf a 2 -f^) + ( ^ ) c 1 ) PQ • 1? 4 6 4 Sin* 

- (A) 15 Sin 4 j. (3 b 4 AoT + ^ b | a a Tv 5 ) (A. 16 ) 


e e dV ± 
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£ c. + a_ b a 4-f4) 1 


e V ± a e e 


e e ^r ± -j 


* l S s a 

+ ( ^ • 24 Sin3 ^ (2 a e |ty 7 b e c e + a e * 3 b 


3 b 
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e e e 


a a 
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^-v n 6 ~ Tvf Cosc< ) 

. a 


+ ^ a e “ b e Cos ^ 3 a e b e c e 


2 ^ o. + a 3 t 3 lie J 


+ a e ■ ^ Tv l 


(A. 17) 


+ 


(terms directly multiplied by c ) 


The partial derivative, given by the Eqs. (A. 8) to 

(A. 17) are calculated at the mean design vector X (that 

is by putting a = a, b = b and c = 0.> 

© 7 0 0 

A. 2 INTERNAL GENEVA MECHANISM: 

C/ A 

In order to calculate ( ) -*■> it is necessary 

i .X. 

to calculate all the quantities on the right hand side 
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at the mean design vector X . From Eq. (2.20), it can 
he seen that 


A. 


b Q Sin c< 

(a 2 + b 2 + 2 a b Cos k )* 5 / 2 
e e e e 


( A . 1 8 ) 


Hence, 


b Sin 


,2 


(A. 19) 


(a 2 + b 2 + 2 ab Cos^ ) 

!1 H! 1 11 

Thus , to evaluate B^ , BJj , B^ , B^ , , A^ and 

a"' at the mean design vector X it is necessary to change 
a to a, b to b and c g to zero. The partial derivatives 

! ?! IT ? ! ?! U ! 

of B 1 , B 1 , B 1 , B 1 , A 1 , Aj , A 1 and Aj with respect to the 
random variable Y ± are given below. 


i a 1 




N V. 

l 


Sin. 


- ( — iy ) . PQ . b SinoC 
2I J ® 


(A. 20) 


where 


Y = 

(a e + 

b 2 + 2 a b 
e e e 

Cos ) 

PQ = 

2 a e 

Ok a 

e 9 h 
?>Y ± + ^ D e 

*< + 


1/2 


(A. 21) 




4* 2 Cos X ( 9 ^ ^ Y + 


'e?-V 


(A. 22) 
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° x 
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+ - rrv Cos<A ) 


f * a e 
+ c e ( TT 7 




+ c \) 

Y 


c e (a e + b e 


Cos ^ ) (-1) PQ 


(A. 23) 
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2*Y^ “ ^ y ^ Cos ^ 4 ' b e Cesd PQ 

, 1 v , ^ a e ,9 . p 

+ ( ? } ( T^ b e + 2 b e TV7 a e } Sin ^ 


+ ( ^> ) a o bf Sin%C (~ 2) PQ 


e e 


(A. 2*0 


B, 


Y 
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* T 1 


t >V 
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*"N £ 

+ (^0 2 Sin-* jja e + b e Cos ^ ) C^- b 0 c e 

^ b e B c e ^ , * a e 
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dV^ e e qY^ e e e e e £ Y . 
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) 1 


(A. 25) 


*" ^ Y ^ 73- V. Sin oC + ( Sin oC ( PQ 
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V — ' 


. t 7* ^ ~r ct £2 u gy — 

rv^e e e a- Y^ 


+ (^) a e b 2 Sin 2 * (- |) PQ 


+ Af) 3 a 2 b^ Sin^ 0 < (- | ) PQ 


+ ^ 3 Sin3 ^ (2 a e |t 7 b e + a e 3 b e |?7 3 


(A. 26) 
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(A. 27) 


- ( ~Y ) 2^® CoscC + (^(l) pq b e cos^ 

+ (^) (->+ a e b 2 Sin 2 ^ + 3 a e b 2 Cos 2 o< )(- |) PQ 
+ (-W) (-^ Sin 2 cx + 3 Cos 2 >( )( 4~w- b 2 

yO C' V * B 


+ 2 a„ b 


zrrpO + 0-rj) 18 a 2 b^ Sin 2 o( Cos U (- |)PQ 


e e ^ v i V e e 
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+ (^-?) 18 Sin^?<, Cos -A (2 a b 3 

y? e e 


+ a e 3 1? SinV a| ^ (- | ) PQ 

+ ( ^7 ) 15 Sln 4 (3 a e It? b e + a e V b l J17 > 

(A. 28 ) 
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The partial derivatives, given by Eqs. (it. 20) to (A. ,29) 

— y~- 

are crlculated at the mean design vector X (that is by 
putting a = a, b = b and c = 0). 
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APPENDIX B 

PROBABILISTIC ANALYSIS OF GEAR TRAIN 
ERRORS 

B.1 INTRODUCTION: 

The primary objective in a precision gear train 
is to provide the functions of motion, power transmission 
and dynamic response with a high degree of accuracy and 
capability. The analysis of gear train performance requires 
identification of error types, sources and their interrela- 
tion to produce a total error. In precision gear trains, 
the errors have to be controlled in an optimal manner. Two 
performance criteria or errors completely define the pre- 
cision of performance of a gear train. These are 

(i) transmission error 
and (ii) backlash. 

(i) Transmission error: 

It is the variation in the transmission ratio of 
a gear pair or train from the ideal nominal value. Although 
in a gear pair each member has a whole number of teeth, 
which fixes the transmission ratio in terms of integral 
numbers of turns, the Instantaneous ratio during one revo- 
lution is likely to fluctuate slightly. Transmission 
error is the measure of instantaneous variation from 
ideal nominal value. This is caused by the net sum of 
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in dividual gear position errors and installation runout 
errors. 


(ii) Backlash : 

It is the total lost motion for a gear pair or 
train. It can also he defined as the amount of motion a 
meshed gear has when its mate is held fixed. It is caused 
by: thinned teeth, enlarged center distance etc. 

Sometimes the net combined effect of transmission 
error and backlash, known as integrated gear train position 
error, is used as a characteristic of performance of the 
gear train. It is the net sum of transmission error and 
backlash and is a measure of total deviation from ideal 
position of a gear pair or train. The various sources that 
cause transmission error and backlash are given below: 


Transmission error sources 


7 

I 


Backlash sources 


1. Position error in the 
individual gears 

(a) Total composite 
error (ICE) 

(b) Supplementary posi- 
tion error (SPE) 

2. Installation errors 

(a) Clearance between 
gear bo:?e and shaft 

(b) Shaft runout at point 
of gear mounting 

(c) Ball bearing radial 
play due to clearan- 

- ces-etc, 


Gear size tolerance and 
and allowance 

Centre distance tolerance 
and allowance 

Eccentricity of fixed bearing 

Ball bearing radial play due 
to clearances 

Shaft runout at point of gear 
mounting 


Thermal dimensional changes 
etc. 
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B*2 TOLEBAUCES FOR C-EAR TRAIN 

Some of the errors that are permitted and are 
given as tolerances in the design of agear pair or train 
are defined below : 

(i) Total composite error (TCE) is defined as the maximum 
variation in centre distance as the gear is rolled, 
intimately meshed with a high quality master gear, on 
a variable centre distance fixture. This is a measure 
of the cumulative effects of eccentricity, tooth to 
tooth variation and profile deviation. Thus TCE cam 
he considered as an overall tolerance and its maximum 
value is specified as x. The permissible values of 
TCE according to AGMA standards are given in 10 . 

The TCE specification can control center distance 
variation, but it is incapable of detecting all angular 
errors. It can reveal tooth-to- tooth profile position 
errors but not accumulated errors. Total angular 
transmission error can only be revealed by checking 
the tooth profile for position error. 

(ii) Gear tooth profile position error or simply the posi- 
tion error is the deviation of any position of a 
gear from the true position it should have. This 

is illustrated in Fig. B.l. This error does not 
become apparent until it is mated in an application 
and until a functional output is observed. However, 
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an individual gear can be evaluated for position 
error by determining the degree of departure of its 
tooth profiles from ideal form and space positions. 

The measurement of the position error can be either 
linear or angular. As a linear measure it is the 
displacement from the ideal position measured along 
the nominal pitch circle. As an angular measure it 
is the subtended angle relative to the center of the 
gear. If the maximum allowable position error is 
specified as + x, then any profile point should not 
deviate from its ideal position in excess of the speci- 
fied amount, x. 

(iii) Centre distance tolerance: 

In practice, the gear center distance will vary 
from the ideal value regardless of the manufacturing 
method used for the gears. The amount of variation 
on tolerance is a function of the fabrication method 
used to establish the gear centers. Tolerance on 
gear center distance is often specified as plus some 
value minus zero, applied to the ideal nominal center 
distance. 

(iv) Gear size tolerance: 

The tooth thickness must be toleranced to per- 
mit realistic fabrication. This tolerance establishes 
an equivalent possible backlash range. The size 
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tolerances are specified as •• + 2 and the values of x for 
various AGHA quality classes are given in 10 

(v) Fixed bearing eccentricity: 

The absolute position of gear centres also 
depends on the tolerances and allox^ances given for the 
installation of bearings. Both ball and sleeve bearings 
that are fixed in a housing by press fit or clamping con- 
tribute to backlash due to the eccentricity that alters 

the effective center distance. The bearing eccentricity 

+ x 

is generally specified as - 0 . 

(vi) Component Shaft run out : 

Shafts, especially if long and slender vill have 
a measurable amount of runout due to boaring (bending). 
(Usually the maximum permissible value (x) of the shaft 
runout is specified). This influences the actual centre 
distance of the gears as shown in Fig. B.2. 

B.3 COMPUTATION OF MEAN AND STANDARD DEVIATIONS OF 
INTEGRATED GEAR TRAIN POSITION ERROR: 

Once the tolerances specified for the gear train 
are known, the mean and standard deviations of the tole- 
ranced quantities can be established (say, by assuming 
three sigma values) and the mean and standard deviations 
■ of transmission error and backlash can be found as indi- 


cated below 
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S t atistical parameters of transmission error t 

(i) F^tablish mean and standard deviation values for each 

transmission error source (e^) as andjT-^j 1 = 1, 2, 


(ii) Find mean and standard deviation of total gear 
transmission error as 


h \ 


O- T 


1 


1/2 r~v~ 


— 1 1/2 


■ 1/2 r^, ,2 o.-,vz 


1 - ( ^ ) 2 


2v w i ) + ( °-i 5 


(B.1) 


Compute jU and a— j using similar formulas 
for the total pinion^ transmission error. 


(iii) Calculate mean and standard deviation of total mesh 
transmission error using the relations 


= 1 <*> 


1/2 r- 2 2 2 2 


1/2 




HT g + M T p 

1 - 


Fg H T p 

for mesh ratio ' 2 * * 5 v 2 
for mesh ratio 2 


1/2 


1/2 


( /% + A, + ^ K + ^J 


“g x "P ' g 
for mesh ratio <C 2 


2 2 1/2 

- (er-m + <r - T ) for mesh ratio ^ 2 

g P 


(iv) Determine the mean and standard deviation of 
entire gear train transmission error as 


(B.2) 



C/Vt =5 (1T > 


i b-i p 


) + J< 


= 2 


^T^m-j 


for stage ratio *C 2 


for stage ratio 2 


v T ; t 


, ** v 

1 - ( T ) 


z< 


‘Mi, 2 V /'’Vm-iPl 
} + 2 ( % ' I 




for stage ratio <C 2 

,V2 

} for stage ratio 2 


(B.3) 


where the subscript (m - i) indicates the ith mesh and 
denotes the stage velocity ratio of ith mesh. 


Statistical parameters of backlash : 


(i) Establish mean and standard deviation values for each 


backlash 


error source (e^) as ^ andcr-^, i = 1 , 2, 


(ii) Find the mean and standard deviation of backlash of 


each gear 


B * °~B 


) and of each pinion 


( Mb , ^ as 

1 p j p 

2 1/2 

^ B g 1 

with similar formulas for Ka ■ „ and^ B 

' P P 


( B .*0 
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(ill) Compute the mean and standard deviation of mesh 
backlash as 


^B^m “ Mb + Mb 

* o ir 

2 2 1/2 

^B^m = ^“B + ^ B ^ 

S P 


(B.5) 


(iv) Determine the moan and standard deviation of back- 
lash of entire gear train as 



Z < 




V 1 


( «Vt 



V i 



(B.6) 


After computing the statistical parameters of transmission 
error and backlash for the entire gear train, the parame- 
ters of the integrated gear train position error j and 
j) can be obtained as folloxrs : 


ytf j = ( / t/ B ) t / 2 + ( // T ) t / Cos 0 

P 2 2 -J 

cr- x = | (cr B ) t / 4- + / Cos 2 0 ! 


(B.7) 


where 0 is the pressure angle 




FIG, B.1 EXAMPLES OF TOOTH PROFILE POSITION ERROR. 
BROKEN PROFILES ARE IDEAL PERFECT-PROFILE 
POSITIONS , SHOWN FOR COMPARISONS WITH THE 
EXAMPLE PROFILES DRAWN SOLID. 

Example \ — A * R ef . Profile ( perfect form and position ) 
B: Spacing error C '.Profile error ( deviation fr6m true 
involute) D . Lead error E 'Radial position error 
(caused by run out ) F-'Tooth thickness error 



FIG. B.2 SHAFT RUNOUT AT GEAR MOUNTING. 



132 


REFERENCES 

1. Hartenberg, R.S. and Denavit, J., Kinematic Synthesis 
of Linkages , McGraw-Hill, New York, 1964- , 

2. Kolhatkar, S.A. and Yajnik, K.S., "The Effect of Play 
in the Joints of a Function Generating Mechanism, " 

Journal of Mechanisms, Vol. 5* No, 4-, 1970, pp. 521-532. 

3. Tuttle, S.B. , "Error Analysis ”, Machine Design, Vol, 32, 
No. 12, I960, pp. 153-158. 

4-. Rothbort, H.A. , Cams-Design, Dynamics and Accuracy , 

Wiley, New York, 1956. 

5. Michalec, G.W., Precision Gearing; Theory and Practice . 
John Wiley & Sons, Inc., New York, 1966 . 

6. Garrett, R.E. and Hall, A. S., "Effects of Tolerances 

and Clearances in Linkage Design," Journal of Engineering 
for Industry, Transactions of ASME, Vol. 91? No. 1, 1969} 

pp. 198 - 202 . 

7. Dhande, S.G. and Chakraborty, J., "Analysis and Synthesis 
of Mechanical Error in Linkages - A Stochastic Approach, " 
ASME Paper No. 72 - Mech. - 36- 

8. Rao, S.S. and Ambekar, A.G., "Mechanical Error Analysis 
of Spherical Function Generating Mechanisms - A 
Probabilistic Approach," ASME Paper No. 74--DET-15. 



133 


9. Dhande, S.G. and Chakraborty, A.J., "Mechanical Error 

Analysis of Cam Follower System - A Stochastic Approach, " 
I. Mech. E., pp. 957-962, 1975. 

10. Kim, H.R. and Newcombe, W.R., "Stochastic Error Analysis 
in Cam Mechanisms", Mechanisms and Machine Theory, 

Vol. 13, No. 6, pp. 631-64-1, 1978. 

11. Rao, S.S., "A Probabilistic Approach to the Design of 
Gear Trains* " Int. <1. Mach. Tool.Des. Res., Vol. 14- , 
pp. 267-278, 1974-. 

12. Dubowsky, S. and F renders tein, F., "Dynamic Analysis 
of Mechanical Systems with Clearances," Journal of 
Engineering for Industry, pp. 305-316,-1971. 

13. Lloyd, D.K. and Lipow, M. , Reli ab ili ty : Management , 
Methods, and Mathematics . Prentice-Hall, Inc., Englewood 
Cliffs, New Jersejr, 1962. 

14- . Rao, S.S., Optimization ; Theory and Applications . 

Wiley Eastern Limited, New Delhi, 1978. 

15. Timoshenko, S. and Young, D.H, and Weaver, W., 

Vibration Problems in Engineering , Fourth Edition, 

John Wiley & Sons, Inc., New York, 1974-. 

16. AGMA 390.02, "Gear Classification Manual," American 
Gear Manufacturers Association, Washington, D.C., 

U.S.A. 



134 


17. Lichwitz, 0., ”1160123111 sms for Intermittent Motion 
Part-I," Machine Design, pp. 134-148, Dec., 1951. 

18. Ray C. Johnson, ir How to Design Geneva Mechanisms 

to Minimise Contact Stress and Torsional Vibration," 
Machine Design, Vol. 28, pp. 107-111, March, 1956. 

19. Amitabha Ghosh and Mallik, A.K., Theory of Mechanisms 
and Machines, Affiliated East-Uest Press Pvt. Ltd., 
New Delhi, 1976. 



